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Abstract 

J/^l' production and collective flow is studied with a coalescence model based on phase space 
distribution of charm quarks from a multi-phase transport model simulation of relativistic heavy 
ion collisions. Both the yield and the flow of J/ip particles are sensitive to charm quark final 
state interactions. As the charm quark rescattering cross section increases from 3 mb to 10 mb, 
J/ip elliptic flow increases faster than corresponding light hadron elliptic flows. The V2{pt) of J/^lJ 
crosses that of D mesons to reach a value that is about the peak value of the D meson flow but 
at a higher pt. As J/V' elliptic flow has only contributions from charm quarks, it complements D 
meson elliptic flow in reflecting charm properties in the Quark-Gluon Plasma. 

PACS numbers: 25.75.-q, 25.75.Ld, 24.10.Lx 
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Experimental data from the Relativistic Heavy Ion Collider (RHIC) at Brookhaven Na- 
tional Lab (BNL) revealed fascinating properties of the Quark-Gluon Plasma (QGP) pro- 
duced in relativistic nucleus-nucleus collisions . The J/ip particle is one of 
the important probes of QGP properties. It was proposed as a signature of QGP forma- 
tion because of its dissociation due to color screening inside the Quark-Gluon Plasma 
Dissociations due to comover scatterings have also been studied to interpret experimental 

QCD) calcu- 
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121. The 



data at SPS energies ja |9|, |10|] . Recently, lattice Quantum Chromo dynamics 
lations show that the J/ip particle can survive the plasma up to about 2Tc 
survival of J/ip was also shown by potential models This leads to new insights into 

the experimental data At RHIC energies, many pairs of charm and anti-charm quarks 
can be produced in a single event. These charm and anti-charm quarks may recombine 
)articles. The recombination can contribute significantly to the final J/ip yield 



into JM particles. 
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2J, |25[ . In this paper, J/ i/j production and flow will be 



studied by a phase-space coalescence model using the charm freeze-out information from a 
multi-phase transport (AMPT) model. In the following, the AMPT model and the coales- 
cence formalism will be reviewed. This is followed by the presentation of results including 
the J/tp yield, (Pt), Pt distributions, and V2{pt)- Finally, a summary will be given with the 
implications of these results on the QGP dynamics. 

The J/ip production in this study is based on the freeze-out phase space information of 



charm quarks from the AMPT mode, 
relativistic heavy ion collisions 



22l 



le AMPT model is a transport model that simulates 
2^ 



0. It uses the HIJING model to provide 
initial conditions. Either mini-jet partons, or partons from string melting will participate 
in the space-time evolution of the system. The parton evolution is carried out by the ZPC 



parton cascade model |3.1|. At parton freeze-out, partons are converted into hadrons using 
either the Lund string fragmentation model or a coordinate-space coalescence model that 

nn 

combines nearest partons into hadrons. Then the ART hadronic transport model [22, 1^ 
is used to evolve the hadronic system. The AMPT model can reasonably describe particle 
distributions at RHIC. It is also successful in showing the importance of partonic evolution 



on elliptic flow and HBT radii 
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In Ref. 133, charm flow in Au+Au collisions 
at ^/sNN = 200 GeV is studied using the AMPT model in the string melting scenario with 
the perturbative method. The initial D mesons follow a parametrization of D meson pt 
distributions from the STAR collaboration j^l and a rapidity plateau between -2 and +2. 
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They are dissociated into charm quarks. Screened Coulomb cross sections are used for the 
rescatterings. Resuhs from the 3mb cross section case are compared to those using a lOmb 
cross section. The D meson elhptic flow and non-photonic electron elliptic flow are found 
to be very sensitive to the charm rescattering cross sections. 

In the following, the charm quark freeze-out information will be used for the study of 



production of the J/ip particle using a phase-space coalescence model 
this model, the J/ip momentum distribution is given by 
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421. In 



J (27r)3x2 /^(^i'Pi)/5(^2,P2)/j/v,(^,9)- (1) 

In the above formula, gj/^ = 1/12 is the degeneracy factor of J/ip production from charm 
and anti-charm quarks, q is the relative momentum, f is the relative position, and R is the 
center-of-mass position of a pair, fc, fs are the freeze-out phase-space distributions of charm 
and anti-charm quarks. /]y^ is the J/ip Wigner function. With a Gaussian spatial wave 
function, 

/j7^(^g)=8expf-4-gV). (2) 



The width a is related to the J/V rms radius by r^^^ = |ct^. Vrms is taken to be 0.5 fm as 
given by the potential model |42| . J/ ip production from two freeze-out distributions will be 
compared. One has 3mb parton rescattering cross sections, and the other has lOmb cross sec- 
tions. Only J/ip production from coalescence is taken into account. There is no production 
and survival from initial nucleon-nucleon collisions. The above coalescence approach can not 
accurately account for binding energy. No feeddowns from higher charmonium resonances 
are taken into account. 

The J/ rapidity density per binary nucleon-nucleon collision as a function of the number 
of participant nucleons is shown in Fig. ^ As expected, the J/ ip yield increases with the 
number of participants. The 3mb case has a larger yield compared to the lOmb case. In 
central collisions, it can be larger by a factor of about 2.4. This is because of charm and 
anti-charm quarks are closer in phase-space with smaller cross section and less rescatterings. 
The solid curves are for the case with 1.4mb for the charm production cross section in 
nucleon-nucleon collisions |4J]. At the moment, the charm production cross section has 
large uncertainty. When the cross section goes down to 0.6mb [4^, the J/ip yield decreases 
to about 20% of the 1.4 mb result. For central collisions, both the 3mb case and the lOmb 



case can give suppression consistent with preliminary PHENIX data 



relative to the 



production in p + p collisions |47[. It is also interesting to see that the centrality dependence 
has almost the same shape as the production from recombination, e.g., from Grandchamp 
and Rapp's calculations {2^. In other words, coalescence and recombination are closely 
related to each other. If the survival of charmonium from initial nucleon-nucleon collisions 
is taken into account, the J/ip results are expected to be similar to those of Grandchamp 
and Rapp. 
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FIG. 1: J/V' yield per binary nucleon-nucleon collision as a function of the number of participant 
nucleons. The data point is from Ref. \^\ - 

The averaged function of the number of participants is shown in Fig. |21 The 10 

mb case has more radial flow and is above the 3mb case. The (pf) increases with centrality 
as more radial flow is generated in the 10 mb case. The 3mb case is slightly different. A 
closer look at the pt distributions reveals that in peripheral collisions, high pt charm quarks 
escape easily and are not affected, while in central collisions, the quenching of high pt charm 
leads to a decrease of {pf) of J/ip particles. The {pf) is about 3 GeV^ in the 10 mb case 
and is about 1.5 GeV^ in the 3 mb. The 10 mb case is comparable to recent calculations 
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from recombinations by Thews et al. 16|, |48| and also comparable to preliminary PHENIX 
central electron arm data 
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FIG. 2: Averaged at mid-rapidity as a function of tiie number of participant nucleons. 

Fig. 121 shows the elliptic flow parameter f 2 as a function of centrality. The shapes of 
the curves are similar to those of charged hadrons. However, different from the charged 
hadron case Q], the elliptic flow is more sensitive to the cross section. This is due to the 
increased sensitivity of radial flow of massive particles relative to light hadrons as seen in 
Fig. El Larger cross section has larger asymptotic elliptic flow at high p^, at the same time 
it has larger {pf) which leads to more weight of high pt flow in the integrated V2. 

More details can be seen by looking at the pt distributions and the pt differential V2 
curves. They are shown in Fig. jHfor minimum-bias Au-|-Au collisions at ^/sNN = 200 GeV. 
In addition to the J/ip results, the charm quark and the D meson (including D* meson) 
results are also shown. Being affected by both the charm and the light quarks, D mesons have 
an invariant pt distribution that has the same concave shape as the charm quark distribution 
with a slightly higher averaged pt- On the other hand, J/ip comes only from charm quarks. 
The Pt distribution has a different, convex, shape, with a more enhanced averaged pt- 
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FIG. 3: Elliptic flow as a function of the number of participant nucleons. 

The elliptic flow results show mass ordering for the low pt region. Charm and D meson 
V2 curves increase with pt, reach a peak, then decrease a little. J ftp V2 in the pt range shown 
here increases, crosses those of charm quarks and D mesons up to a value that is comparable 
to the peak value of charm and D meson elliptic flow. This behavior is different from some 

n 

previous studies in which J/ip V2 is consistently larger than D meson V2 |42|. The crossing 
of the J/ip and D meson V2 curves reflects the distinct freeze-out phase space distributions 
from the AMPT model. 

In summary, J/ ip production from the coalescence of charm and anti-charm quarks reflects 
charm interactions in the QGP. The J/ ip yield is very sensitive to the charm production cross 
section in p+p collisions. Coalescence yield has a centrality dependence that is similar to that 
of recombination models. Both the 3mb and lOmb cases can lead to suppression on the same 
level as recombination models and consistent with preliminary PHENIX results. The lOmb 
case gives (p'f) around 3 GeV^ and is comparable to recombination results and preliminary 
PHENIX data, while the 3mb results are much lower than the lOmb case. The elliptic flow 
follows the same centrality dependence as charged hadrons, but is more sensitive to the cross 
section partly due to the enhanced sensitivity of pt distribution because of the large mass 
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FIG. 4: Invariant transverse momentum distributions and pt differential elliptic flows for charm 
quarks, D mesons, and J/V' particles. 

of the J/ ip particle. The J/ ip pt distribution is convex and different from that of D mesons 
because of the combination of charm and anti-charm quarks. The J/ip elliptic flow crosses 
that of D mesons and reaches about the peak value of the D meson elliptic flow. Hence, the 
combination of open and hidden charm measurements will provide a more complete picture 
of the evolution of the charm quarks in the strongly interacting Quark-Gluon Plasma. 
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